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ABSTRACT 

In this paper we present the coordinates of 67 55' x 55' patches of sky which have the rare 
combination of both high stellar surface density 0.5 arcmin -2 with 13 < R < 16.5 mag) and 
low extinction (E(B — V) ^ 0.1). These fields are ideal for adaptive-optics based follow-up of 
extragalactic targets. One region of sky, situated near Baade's Window, contains most of the 
patches we have identified. Our optimal field, centered at RA: 7 h 24 m 3 s , Dec: — 1°27'15", has an 
additional advantage of being accessible from both hemispheres. We propose a figure of merit for 
quantifying real-world adaptive optics performance, and use this to analyze the performance of 
multi-conjugate adaptive optics in these fields. We also compare our results to those that would 
be obtained in existing deep fields. In some cases adaptive optics observations undertaken in the 
fields given in this paper would be orders of magnitude more efficient than equivalent observations 
undertaken in existing deep fields. 



Subject headings: Astronomical Techniques, Astronomical Instrumentation, Astrophysical Data, Galaxies 
1. Introduction 



1 Department of Astronomy & Astrophysics, University 
of Toronto, 50 St. George Street, Toronto, ON, M5S 3H4 
Canada 

2 Centre for Astrophysics and Supcrcomputing, Swin- 
burne University of Technology, 1 Alfred St, Hawthorn, 
Victoria 3122, Australia 

3 Research School of Astronomy and Astrophysics, In- 
stitute of Advanced Studies, The Australian National Uni- 
versity, Canberra A.C.T., Australia 

4 Gemini Observatory, Southern Operations Center, c/o 
AURA, Casilla 603,La Serena, Chile 

5 Observatories of the Carnegie Institution of Washing- 
ton, 813 Santa Barbara Street, Pasadena, CA 91101 

6 Sterrewacht Leiden, RO. Box 9513 NL-2300 RA Lei- 
den, The Netherlands 

7 Centro de AstroHsica, Universidade do Porto, Rua das 
Estrelas, 4150-762 Porto , Portugal 

8 Canada- France-Hawaii Telescope Corporation, 65-1238 
Mamalahoa Highway, Kamuela, Hawaii 96743, USA 

9 Institut d'Astrophysique de Paris, UMR7095 CNRS, 
Universite Pierre et Marie Curie, 98 bis Boulevard Arago, 
75014 Paris, France 

10 The United States Naval Observatory, 3450 Mas- 



Our understanding of the high-redshift uni- 
verse has been revolutionized by deep fields, sev- 
eral of which have been extensively surveyed at 
all accessible wavelengths. Figure Q] shows an 
up-to-date summary of the locations of all ex- 
isting deep fields (red circles). These fields have 



evolution out to very high redshifts (Cowie et al. 




1995 




Yee et al.ll2000: 


Labbe et al.ll2003: Bell et al. 




2004 




Eeami et al.ll2004t Ivan Dokkum et al.l 20041 


Arnouts et al.ll2007; Davis et al.l2007llScoville et al. 


2007; Bouwens et al. 


2010; Ono et al.l 2010h. Be- 



cause galaxies at such high redshifts are typi- 
cally < 1 arcscc in size, kinematical investiga- 
tions of galaxies in these fiel ds require adap- 



tive optics (AO) spectroscop y ( Law et al. .2009; 



2PJ 

Forster-Schreiber et al.l 12009). The promise of 
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Fig. 1. — The galactocentric coordinates of existing deep fields (red circles) and the locations of the fields 
better suited for AO observations presented in Table 1 of this paper (blue squares, see Section [3] for details). 
The location of the fields has been overplotted on the dust emission map from the Schlegel, Finkbeiner & 
Davis (1998) study. Labeled existing deep fields are indicated with filled circles. The green square denotes a 
1 square degree region within the CFHTLS W2 field (Section l5.2|) . The field labeled 'AODF' is our suggested 
optimal field whose properties are studied in detail in Sections |4] and [5] 



such observations has been held out as an excitin g 
next step for over a decade fe.g.. lEllis et al.lll997[ ). 
Unfortunately, it is now clear that only very lim- 
ited AO observations are going to be undertaken 
in any existing deep fields. 

Coupling integral-field spectroscopy to AO is 
crucial for understanding the formation of massive 
galaxies, particularly disks, since at high-redshift 
it has proven difficult for slit spectroscopy to re- 
liably identify kinematic disks as kinematic and 
morphological ax es are not necessarily correlated 



rotati onal gradients ab ou t a preferred kinematic 



(|Erb et all 12003 ) . Even at intermediate redshift 



(z ~ 0.6), it has been demonstrated that galax- 
ies are already kinematically complex and that 
3D integral field spectroscopy (IFS) is essential 
to p hysical understanding and kinema tic model- 
ing (|Rix et all Il997t iFlores et all l2006h . At the 
highest redshifts AO IFS observations by some 
groups have given different results compared to 
non-AO observations of other groups. For ex- 
ample Laser Guide S tar (LGS) AO o bservations 
with kpc resolution ( Law et al. 20091 ) show that 
z = 2 — 3 Lyman-break selected galaxies have high 
intrinsic velocity dispersions and no significant 



axis 



(|Law et al.1 120091 ). iForster-Schreiber et al 
( 20091 ) found similarly high velocity dispersions 
but a much greater incidence of disk rotation in a 
predominantly non-AO dataset of predominantly 
near-IR selected galaxies. It remains inconclusive 
whether such differences arise from a difference in 
sample (massive vs low mass galaxies) or the fact 
that the non-AO data has seven times poorer res- 
olution on average in natural seeing. This is an 
important question: physical differences in kine- 
matics at high-redshift may diagnose the preva- 
lence of fast gas accr etion along cold flows in the 



early Universe (e.g., iBournaud et al.l 120071 ). but 



they may also arise from sample select ion effects or 



observ ational limitations (for example lGreen et al 



(2010) suggest it is simply the high-star-formation 
rates which drive the large velocity dispersions). 
Most IFS observations at high-redshift are still 
done without AO due to the technical difficulties 
of AO and also to the practical difficulties of find- 
ing enough targets near sufficiently bright stars in 
existing deep field samples. 

The next stage in the development of this field 
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is to complement kinematical studies by probing 
chemical abundance gradients at a sub-kpc scale 
in star-forming galaxies, and to extend existing 
kinematical investigations to encompass more rep- 
resentative galaxies. This requires AO systems to 
be operating more efficiently (i.e. without per- 
formance limitation imposed by natural guide- 
star availability) and, ultimatelely, to multiplex 
if truly large samples are to be obtained. A step 
in this direction is already being taken with the 
MASSIV survey on the VLT, which targets star- 
forming galaxies in the redshift range 1 < z < 2 
with SINFONI (|Epinat et al.l 120091 : iQuevrel et al ' 



(|2007l ) a re accessible to the V LT laser guide star 



system ( Bonaccini et al. 19991 ): the loss of 75% 



2009). The targets are more representative than 
those being probed by SINS, with median stellar 
masses of ~ 1O 1O M0 and median star formation 
rates of ~ 10M Q yr _1 (the corresponding values 
for SINS are - 10 la5 M Q and - SOMoyr" 1 , re- 
spectively). However, most of the SINFONI data 
acquired during the MASSIV survey are seeing- 
limited leading to a final median spatial resolution 
of - 0.6-0.7" with only 25% of the MASSIV sam- 
ple presently being observed with adaptive optics. 
Of these AO targets only a few are being acquired 
with the smallest pixel size (0.05"). The main rea- 
sons are (i) the limitations in the availability of 
natural guide stars which precludes usefully ob- 
serving at finer available pixel scale, and (ii) the 
difficulty to reach the depth required to probe the 
low-surface brightness component of galaxies in a 
reasonable exposure time with the smallest pixel 
size. This latter point leads to expectations of con- 
siderable progress in this subject with the advent 
of 30m-class Extremely Large Telescopes (ELTs). 

A basic problem with undertaking AO in ex- 
isting deep fields, even with laser guide stars, is 
that one still needs at least one reasonably proxi- 
mate natural guide star to sup ply the information 
neede d for tip-tilt correction (jRigaut fc Gendro"nl 
1992). In contrast, two of the main selection cri- 
teria when identifying deep fields have been that 
they contain as few bright stars as possible to 
avoid light scattering contamination and satura- 
tion in long exposures, and that th ey lie in re- 



gions of low Galactic extinction (e.g.. lAlcala et al 



20041 ). Thus all existing deep fields are near the 
Galactic poles, where the density of suitable nat- 
ural g uide stars is near a minimum. For ex- 
ample, iDavies et al.l ([20081) report tha t only 1% 
of the Lyman break galaxy sample of iMannucci 



of the targets is due to the absence of suitably 
close natural guiding stars (NGS), while additional 
25% are lost after suitable color cuts and elimi- 
nation of systems at redshifts obscured by strong 
OH features. The situation is similar with Gem- 
i ni, whose AO system h as similar sky coverage 
([Ellerbroek fc Tylerlll998l ) . Even with the upcom- 
ing Gemini Multi-Conjugate AO system (MCAO), 
the -ff-band sky covera ge at the galactic poles 
will only be around 15% ( Rigaut et al.llio OQ) . and 
large benefits for MCAO emerge from having more 
than the minimum number of natural guide stars. 
This is because the geometry of the guide stars 
on the sky impacts the un iformity of Strehl ra- 
tio ([Flicker fc Rigautll200l[) . 

The issue of guide star rarity in deep fields be- 
comes prominent in cases where target source den- 
sity is low. This is often the case for extragalactic 
programs which focus on unusual objects. For ex- 
ample, many of the key projects d escribed in the 



JWST Design Reference Mission ([Gardner et al 



2006) rely on cither extreme depth or serendipi- 
tous lines-of-sight. If such JWST observations are 
to be synergistic with ground-based AO follow- 
up, in particular with next generation telescopes 
like TMT or E-ELT, they cannot be undertaken 
efficiently in any existing deep field. It would 
be disappointing indeed if only 1%-10% of rare 
targets imaged with JWST in a deep field could 
be followed-up with a ground-based integral field 
units (IFU). It is becoming clear that existing 
and planned AO systems are set to enable trans- 
formative high-redshift science, but they will do 
so only in the regions of the sky in which they 
are effective. It is arguable that no existing deep 
field is suitable for efficient extragalactic AO work 
(though of course the cost of obtaining ancillary 
data equivalent to that already obtained in exist- 
ing deep fields may overwhelm the gains obtained 
from high-efficiency adaptive optics). 

In this paper we report on the results we 
have obtained in searching for those fields on the 
sky most suitable for high-efficiency extragalactic 
adaptive optics observations. In Section [2] we de- 
scribe the important characteristics of deep fields 
in the context of adaptive optics observations, 
such as the acceptable level of dust extinction, 
field size, and magnitude range of natural guide 
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stars. In Section [3] we describe our attempts to 
identify the most suitable areas on the sky for un- 
dertaking extragalactic AO work, which is based 
on the strategy of combining information from 
all-sky stellar density and extinction maps. Our 
preferred 'AO-friendly' field and its first imaging 
results are described in Section 21 In the following 
Section [5] we define a figure of merit for adaptive 
optics and use this to compare the efficiency of AO 
observing in the proposed fields relative to the ef- 
ficiency in representative deep fields. Our conclu- 
sions are summarized in Section |6l All magnitudes 
in this paper are based on the Vega system. 

2. Desirable characteristics of extragalac- 
tic fields optimized for adaptive optics 

In this section we consider the desirable charac- 
teristics of extragalactic fields optimized for adap- 
tive optics. The relevant considerations include 
the maximum acceptable level of extinction from 
the intergalactic medium, the minimum useful 
area on the sky, and number density and mag- 
nitude range of available natural guide stars. We 
will consider each of these factors in turn, and dis- 
cuss the importance of each of these factors using 
very general principles, in order to look for consid- 
erations that will remain relevant for future AO 
systems. 

2.1. Extinction 

Although most adaptive optics is undertaken in 
the near-infrared where extinction is lower than at 
visible wavelengths, it is clear that for any num- 
ber of reasons, including reliability of photometric 
redshifts and 'future-proofing' the fields so as to 
make them useful when AO work moves to shorter 
wavelengths, that the ideal fields will lie in re- 
gions of low Galactic extinction. As any glance 
at the night sky will attest, patchy extinction can 
be rather high in regions with high star counts. It 
is therefore important to define an upper limit to 
the acceptable extinction in order to exclude un- 
suitable fields. A value of E(B — V) ~ 0.15 mag 
i s a g ood s tarting point, becau se iFukugita et al 



(|2004) and lYasuda et all (|2007l ) show that galac- 
tic extinction estimates become fairly unreliable 
in regions with E(B — V) > 0.15 mag. To err 
on the conservative side, in this paper we will use 
an upper limit of E(B — V) = 0.1 mag on the 



mean extinction as a constraint when exploring 
star count surface density maps for suitable fields. 
We note that E(B — V) — 0.1 mag corresponds 
to A v — R v x E(B - V) ~ 0.3 mag at visible 
wavelengths, and that this is a factor of three to 
ten times higher than the corresponding extinction 
in near-infrared (NIR) passbands used by current 
AO systems. 

2.2. Field size 

The next factor to consider is the required size 
of the field. For extragalactic fields, the area of 
the field is driven by a desire to minimize the im- 
pact of cosmic variance, because scale-dependent 
inhomogeneity is often the dominant source of er- 
ror in measurements derived from galaxy popula- 
tions within a survey volume. The survey volume 
naturally depends on the area on the sky and the 
chosen redshift range, but for concreteness we will 
assume that most extragalactic work will explore 
a range of redshifts from z = to z = 4, which 
encompasses most of the star-formation history of 
the Universe. For such surveys, areas on the order 
of a square degree are needed in order to maintain 
fractional errors on number counts near the 10% 
level, and to probe a wide range of cosmic struc- 
tures. This is fairly easy to demonstrate using on- 
line tools su ch as the Cosmic Varianc e Calculatoi0 
described in lTrenti &: Stiavelli but an even 



simpler way to show t his is to use the analytic ex - 
pressions provided by iDriver fc Robothanl (l2010h 



to estimate and compare cosmic variance for differ- 
ent field sizes. These authors employed counts of 
galaxies near the characteristic break in the lumi- 
nosity function (M*) in the Sloan Digital Sky Sur- 



vey D ata Release 7 (SDSS DR7, lAbazajian et al 



2009) to derive an empirical expression connecting 
cosmic variance and survey volume. Assuming a 
single sight-line and a rectangular geometry, the 
fractional error in the counts of M* galaxies is 
given by: 



(1- 0.03(^04/5) -1) 

x (219.7 - 52.4(log 10 [A ■ B ■ 291.0]) 

+3.21(log 10 [A • B • 291.0] 2 )) 



^(7/291.0. 



(1) 



http: //casa. Colorado . edu/~trenti/CosmicVariance . 
html 
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Fig. 2. — Cosmic variance, quantified using Eq. [T] as a function of redshift for four fields covering different 
areas on the sky. Redshifts of presented points correspond to median values for redshift ranges indicated 
with doted lines. The effect of small-scale inhomogeneity on the field size we propose (~ f square degree, 
denoted as AODF) is comparable to the COSMOS field cosmic variance, and much less prominent than in 
the other two (smaller) fields, GEMS and GOODS. 



where A, B, and C are the median redshift trans- 
verse lengths and the radial depth of the sur- 
vey, respectively, expressed in units of h^ j Mpc. 
(Note that the derived cosmic variance is for 
M* ± 1 mag population only and will take 
higher values for more massiv e halos, see e.g., 
Moustakas fc Somervillle 2002). Results com- 
puted using this equation are presented in Fig- 
ure^ which shows the calculated cosmic variance 
for a number of surveys, and compares these with 
our proposed field size of around one square degree 
(actually 55' x 55', for technical reasons described 
below). 

Figure [5] shows that the calculated cosmic vari- 
ance for our proposed field size results in frac- 
tional counting errors of around fO — 15% (per 
unit redshift interval) for counts of M* galax- 
ies at redshifts between z = 1 and z = 4. 
This is onl y slightly higher tha n for the COS- 



MOS field (|Scoville et al.l 12007ft . but quite sig- 

smaller volume sur- 
and 



nificantly better than for 
vey fields, such as GEMS (IRix et al. 
GOODS (|Dickinson et al.l 120031 ). On this ba- 



2004) 



sis alone we would argue that something around 
square degree is probably the right minimum size 



for a contiguous area survey field intended to al- 
low a broad range of investigations using adaptive 
optics, although another important factor is that 
a survey of this size will contain many thousands 
of strong line emitting objects, which are obvious 
targets for present-generation AO systems. 

We have computed the surface density of 
strong H Q line emitters (which we define to 
be Fjj a > 10~ 16 erg cm 2 s -1 , corresponding 
to the flux density of brig ht line emitters in 



Fdrster-Schreiber et al. 2009; Law et al. 2009) on 



the basis of direct measurement ( Villar et ail 20081: 
Shim et al.1 [2009) as well a s using indirect esti- 



mates scaled from UV flux (Bouwcns et al 



and measurements of [Oil] ([Cooper et al 



2009) 
2008). 



By incorporating all the available information we 
estimate this value to be 2-5 H a line emitters with 
flux > 10~ 16 erg cm 2 s _1 A -1 per square arcmin 
at 1 < z < 1.5, declining to 1-2 per square arcmin 
in the redshift interval 2 < z < 2.5. The deep 
fields proposed in this paper will thus have around 
10,000 suitable targets for AO-based follow-up. A 
significant fraction of these will be lost for various 
reason s (e.g., if H a lie s on an airglow emission 



line, iDavies et al.l 120081 ). and a small number of 
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remaining objects will still lack suitable guiding 
stars (see Section [5J. However, thousands of AO- 
accessible targets will remain, presenting a mul- 
tiple order-of-magnitude change from the current 
situation. 

2.3. Guide star limitations 

We now explore the brightness of natural guide 
stars needed for effective use of adaptive optics. 
Our focus will be on the the following three classes 
of AO systems: 

• Case 1: Laser-assisted adaptive optics sys- 
tems on 8m-class telescopes, for which nat- 
ural guide stars are needed to supply tip- 
tilt corrections. Such systems will define the 
state of the art for the next few years. 

• Case 2: Ground-layer adaptive optics sys- 
tems for 4m-class telescopes. Such sys- 
tems are now being proposed as a means 
of revitali zing 4m-class fa cilities (e.g CFHT 



'IMAKA. lLai et al.1 12008). These facilities 



will also require natural guide stars for tip- 
tilt correction. 

• Case 3: AO systems on 30m-class telescopes, 
some designs for which rely on AO for rou- 
tine operation. In this case we mainly seek 
fields with an abundance of natural guide 
stars bright enough to feed laser-unassisted 
AO systems. Laser beacons may not be 
available at all times, and the existence of 
extragalactic fields in which they are not es- 
sential may be extremely attractive for tele- 
scopes that heavily emphasize AO. 

We will begin by first outlining the general 
problem before focusing on the parameter space 
appropriate to the specific cases above. As will be 
shown below, in practise it is Cases 2 and 3 that 
drive our chosen magnitude limits. 

In order to function an AO system needs to 
capture photons from a star, compute a correc- 
tion, and apply this correction to an optical sur- 
face. The frequency over which an AO system 
must operate is set by the velocity of the atmo- 
sphere and the atmospheric coherence length. The 
coherence length is the length scale over which the 
index of refraction of the atmosphere is effectively 
constant, and is typically around 10 cm at a good 



site. Wind speeds in the upper atmosphere are 
around 20 m/s, so it typically takes around 0.005s 
for a p atch of atmosphere to m ove a coherence 
length ( Rogemann fc Welshlll996 ). The minimum 
frequency of an AO system is therefore around 
200 Hz, although in reality one would want to both 
Nyquist sample the signal and allow time for ac- 
tuator lag in applying a correction, so a realistic 
minimum is around 1 kHz. 

How many photons from a natural guide star 
are needed in this time depends on the specific 
type of correction, but we can bracket our analy- 
sis by considering two extremes: (i) tip-tilt correc- 
tion, for which relatively faint stars suffice, and (ii) 
full correction to obtain diffraction-limited perfor- 
mance, for which bright stars are needed. 

A zeroth magnitude star has a i?-band flux of 
3080 Jy at the top of the atmospher cl, correspond - 
ing to 2.02 x 10 11 photon/m 2 /s (jBesselll [l979h . 
Thus an 8m telescope captures ~ 600 i?-band pho- 
tons from an 18th magnitude star in one millisec- 
ond. In the foreseeable future no AO system will 
have a quantum efficiency approaching unity, but 
even with an end-to-end efficiency of 20% over 100 
photons will remain, which is ample for obtaining 
a reasonable centroid. Thus, at least in principle, 
an AO system on an 8m telescope can use R = 18 
mag stars for tip-tilt corrections. Since the num- 
ber of photons from a star imaged with a 30m 
telescope is about fourteen times greater than for 
an 8m telescope, a 30m telescope can do tip-tilt 
corrections on guide stars down to around R = 21 
mag. On the other hand, a 4m telescope needs 
stars of about R = 16.5 mag for tip-tilt correc- 
tions. We emphasize that these numbers are all 
for rather idealized AO systems. For example, in 
the real-world situation of the Gemini Altair AO 
system, tip-tilt reference stars of around R ~ 16.5 
mag (over a magnitude brighter than the some- 
what ideal case discussed above) are found to be 
highly desirable for high-performance AO opera- 
tion. 

Much brighter natural guide stars are needed 
for use with natural guide AO star systems that at- 
tempt to achieve diffraction-limited performance. 
In this case the size of relevance is not the full aper- 



2 For concreteness we consider the brightness of guide stars 
at visible wavelengths, though the argument can be gener- 
alized to stars at arbitrary wavelength. 
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ture of the telescope, but rather the sub-aperture 
defined by the coherence length of the atmosphere 
which in turn drives the number of needed actu- 
ators. An R — 13 mag star supplies ~ 10 pho- 
tons in 1 ms to a 10 cm diameter sub-aperture. 
The number of photons per bin needed to reliably 
compute a wavefront depends critically on factors 
such as the read noise of the detector, but ten 
photons per coherence-length sized patch on the 
pupil is a reasonable lower limit. Note that in 
the case of diffraction-limited AO (and unlike the 
case with tip-tilt correction), having a larger tele- 
scope does not gain one a fainter magnitude limit 
for natural guide stars, and in fact AO becomes 
harder because the system requires more actua- 
tors. We also note that the presence of bright stars 
(R < 14 mag) in the field with AO correction can 
potentially cause problem for infrared (IR) detec- 
tors by leaving long-lasting (up to an hour) resid- 
ual flux. Although this would affect imaging in 
case when AO correction is applied across a wide 
field of view, the main motives for developing an 
AO-friendly deep field (see Sections Q] and [5]) are 
high resolution IFU or multi-object spectroscopic 
surveys, that would not be influenced. 

On the basis of the considerations just given, 
our search for locations on the sky suitable for ex- 
tragalactic adaptive optics focuses on stars in the 
magnitude range 13 < R < 16.5 mag. The bright 
end is set by the apparent magnitude of stars 
needed to supply guide stars for natural guide-star 
AO systems (independent of telescope aperture), 
which is essentially Case 3 above. The faint end is 
set by the apparent magnitude of stars needed to 
supply tip-tilt reference stars for real- world opera- 
tion of existing AO systems on 8m-class telescopes 
(Case 1) and for ideal-case laser-based ground- 
layer AO with 4m-class telescopes (Case 2 above). 

3. Identification of the suitable fields for 
Adaptive Optics 

In order to find the regions on the sky with the 
properties we have just described, we rely on full 
sky reprocessed composites o f the COBE/DIRBE 
and IRAC/ISSA dust maps (jSchlegel et al.lll998h 
and the UCAC2 astrometric catalog of ~ 5 x 
10 7 stars with d eclination in the [— 90 deg, +(40 — 
52)deg] range (jZacharias et al.l 120041 ) . We con- 
structed a full sky map of star count surface 



density for UCAC2 stars in the range 13 — 
16.5 mag using the HE A LP IX data analysis pack- 
age (jGorski et al.1 120051 ) that performs pixeliza- 
tion of the sphere with equal area pixels. Two 
maps have been produced: one with the reso- 
lution of 6.'871 to match the resolution of the 
available HEALPIX map of Galactic reddening 
E(B — V) and the other (see Figure [3]) with the 
coarser sampling of 55' (the HEALPIX resolution 
that is closest to the 1° x 1° field size, see Sec- 
tion ^. 2p . The resolution of the existing E(B — V) 
map was degraded to match the 55' resolution of 
the star counts surface density map by taking the 
average extinction value for each cell. The coarse 
resolution extinction and star count maps are both 
shown as panels at the top-left and top-right of 
Fig. [3] Note that the UCAC2 catalog has a gap 
in coverage at high declination (shown in gray in 
the figure), but any AO-optimized fields which 
might exist at these very high declinations would 
be generally unsuitable anyway. Any such fields 
would be inaccessible from Chile and be at quite 
high airmass most of the time for major north- 
ern hemisphere observatories (including those on 
Mauna Kea). 

Before proceeding with a detailed analysis, it 
is instructive to note that many positions in the 
sky likely to be suitable for our purposes can be 
identified easily by simply looking for maxima in 
a map obtained by multiplying the stellar density 
map by the inverse of the extinction map. This is 
shown as the large bottom panel in Figure[3l Local 
maxima in this map do not necessarily define re- 
gions suitable for AO, because some local maxima 
correspond to regions with low star counts but ex- 
tremely low extinction. However, this figure acts 
as a natural starting point for the next step in our 
analysis. 

Having identified candidate fields using the 
analysis just described, we then looked at all 
the candidate fields individually to try to bet- 
ter understand their characteristics. To be ex- 
plicit, we first identified all HEALPIX cells whose 
13 < R < 16.5 mag stellar density S sc was 
E sc > 0.5 arcmin~ 2 and whose extinction was 
E(B — V) ^ 0.1. We found 442 one square degree 
cells met these criteria, and these were then ex- 
amined further. The distribution of stellar density 
and extinction for these cells is shown in Fig. |4j 
color-coded by right ascension. In order to cull 
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Fig. 3. — (Top left:) An all sky map of extinction, scaled logarithmically. The solid line grid corresponds to 
the celestial coordinate system with RA in degrees increasing to the left. Zero degrees lies at the center of 
the figure. A Galactic coordinate system is over-plotted with dashed lines. (Top right:) The corresponding 
map of star count surface density for stars in the 13 — 16.5 magnitude range. The region shown in gray 
corresponds to a high declination gap in coverage in the UCAC2 stellar catalog. As noted in the text, 
any AO-friendly fields which might exist at very high declination would be unsuitable for other reasons. 
(Bottom:) A map constructed by multiplying the map at the top left by the inverse of the map at the top 
right. Maxima in this figure correspond to potentially interesting locations for undertaking extragalactic 
adaptive optics observations. Red circles present the positions of 67 fields well-suited to extragalactic AO. 
See text for details. 
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Fig. 4. — Extinction E(B — V) as a function of the star counts surface density E sc for 442 55' x 55' fields with 
S sc > 0.5 arcmin -2 and E(B — V) ^ 0.1. The fields are color-coded based on their equatorial coordinates. 
The dashed line encloses 67 fields with S sc ^ 0.65 arcmin^ 2 and 0.05 < E(B ~ V)[mag] < 0.087. The fields 
flagged with open circles have the highest star counts surface density or the lowest mean extinction or its 
standard deviation. Colored arrows point at the representative fields for each of the three sightlines (see 
Appendix lAl for details). The proposed 'optimal' field described in Section [4] is labeled 'AODF' and flagged 
with an open box. 
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Fig. 5. — The three higher-order moments of the extinction and star count surface density distributions as 
functions of the mean star count surface density for 67 fields from Table [1] The fields are color-coded based 
on their equatorial coordinates, as given in Figure |4j The fields flagged with open circles or with colored 
arrows correspond to the flagged fields in Figure |U Our optimal AO-friendly field is labeled as in Figure |U 
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these fields down to a more manageable number, 
we then restricted the sample further to include 
only those fields whose extinction is lower than 
the average extinction, and whose stellar density 
is higher than the average density. This trans- 
lates into selecting fields with a stellar density 
greater than 0.65 arcmin and reddening less 
than ~ 0.087 mag. This final cut corresponds to 
selecting fields inside the dashed region shown in 
Fig. SI and brings the total number of fields down 
to 67 from 442. The positions of these fields are 
shown as solid dots in Figure [3] and blue squares 
in Figure [TJ 

The locations and properties of the 67 fields are 
given in Table [T] Variances and higher-order mo- 
ments of stellar counts and extinction within each 
field are also tabulated. These moments are based 
on computations within sub-cells with widths of 
6.'871. The fields flagged with open circles in Fig. [4] 
have the highest star count surface density or the 
lowest average value or standard deviation of the 
reddening coefficient, and rows corresponding to 
these fields are italicized in the table. Colored 
arrows point at some more representative fields 
for the three sightlines (see Appendix [A] for de- 
tails). A graphical summary of all higher-order 
statistics is presented in Figure [5j from which it 
can be seen that there is a substantial variation 
in the distribution of both stars and extinction 
throughout the fields we have identified. While 
all the tabulated fields should be quite good for 
extragalactic adaptive optics work, the best field 
for a given purpose will depend upon the specific 
application (e.g. upon whether it is more impor- 
tant to maximize uniformity of star counts, or to 
minimize absolute extinction, or to minimize vari- 
ation in extinction). The various trade-offs that 
need to be balanced in order to choose the best 
field for a given set of requirements are explored 
in Appendix [Al 

4. CFHT imaging results for a proposed 
'optimal' field 

All fields labeled in Figure [1] and in the lower 
panel of Figure [3] have exceptional characteristics 
in terms of stellar surface density and/or extinc- 
tiorE However, a practical factor that has not 



yet been considered is the position of the field 
for accessibility with a broad range of telescopes. 
An ideal field lies near zero degrees declination 
in order to be reachable from both hemispheres. 
With this consideration in mind, the most inter- 
esting region for further analysis proves to be the 
rather large 30° x 60° region centered at RA = 8 h , 
DEC = 0° (see the bottom map in Figured]). To 
identify the best one square degree patch within 
this region, we tessellated the region into 55' x 55' 
cells and explored the distribution of star counts 
and E(B — V) on a cell-by-cell basis. After identi- 
fying a handful of promising cells, we then looked 
for both low absolute extinction and uniformity 
in extinction within individual cells. (Uniformity 
in extinction is desirable for accurate photometric 
redshifts). Figure [6] shows the intra-cell stellar sur- 
face density and the extinction map for the best 
55' x 55' cell, which we refer to as the 'Adaptive 
Optics Deep Field' (AODF). The field lies at RA: 
7 h 24 m 3 s , Dec: -1°27'15" and is labeled with ID 
8328 in Table [1] It has a stellar surface density of 
more than two stars per 2' x 2' region over > 99% 
of its area. In addition, for > 85% of the field pre- 
sented in Figure |6] extinction is E(B — V) ^ 0.1 
(Ay < 0.3; NIR extinction at AO wavelengths 
will be far lower, see Section [2]). Another impor- 
tant practical consideration is that the number of 
very bright stars (which scatter and raise the sky 
background) in this ~ 1 square degree field is low: 
there are only a handful of stars brighter than 11th 
magnitude in the field. 

Since this 'AODF' field seems highly interest- 
ing for future follow-up, Director's Discretionary 
time on the Canada-France-Hawaii Telescope was 
used to explore its properties further (and to act 
as a sanity check on the analysis presented in this 
paper). A 10 min snapshot in g' and z' bands was 
acquired using CFHT's MegaCam in March 2010 
in order to evaluate the distribution and color of 



the b rightest stars in the AODF (see lBertin et al 



120021 for the description of the TERAPIX software 
modules designed for processing MegaCam data). 
The upper panel of Figure [7] presents the ~ 1° x 1° 



3 It is interesting to speculate on why these fields exist. 
While many of them lie in the vicinity of Badde's window, 



characterized by very low dust content, a smaller group of 
fields is found near the Galactic anti-center. One possibil- 
ity, suggested to us by Sidney van den Bergh, is that these 
fields have high counts but low extinction because metallic- 
ity decreases with increasing Galactoccntric distance. As a 
result the Galactic anti-center has a low dust to gas ratio. 
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Fig. 6. — Extinction and star surface density maps centered on the position of our chosen field. Each pixel 
is 6.'871 x 6.'871. There are > 2200 13 — 16.5 mag stars suitable for AO guiding in the field corresponding to 
~ 3300 star- forming galaxy candidates with the number surface density of 3 arcmin -2 . The field is centered 
at RA = 7 h 24 m 2.67 s , DEC = -1°27'14.44". Solid line grid corresponds to the celestial coordinate system 
with RA[°] = 360° - a[°] 



field of view (made out of 5 dithered exposures), 
and it is apparent here that bright stars are indeed 
sparse, with the most prominent one (V=9.9 mag) 
in the center of the field. The main drawback of 
such a star is not the vertical blooming which af- 
fects a small fraction of the imaging area but the 
halos due to internal reflections in the MegaPrime 
optics: such a halo increases locally the sky back- 
ground and limit detectivity. The upper panel of 
Figure [7] shows that there are four stars that cause 
a potential problem. Each halo covers a disk of 
3' in radius, leading to ~ 120 square arcmin for 
the whole field. When compared to the MegaCam 
field of view of ~ 1 square degree, those four halos 
produce a negligible loss of less than 4%. (We note 
that brighter galaxies can still be extracted from 
these areas). 

We also compared AODF galaxy counts in g' 
and z' bands with the expected number based on 
the CFHT Legacy Survey (CFHTLS) Deep data 
(details on galaxy/star separation method used by 
the T ERAPIX pipeline are given in lCoupon et al. 
2009). The resulting depth found in both g' and 
z' bands follow the expectations (within the range 



of error), and the recovered galaxy number surface 
density closely tracks the distribution of CFHTLS 
Deep objects, as shown in the lower panel of Fig- 
ure [Jj The slight excess at the bright end is the 
combined effect of shot noise and cosmic variance 
(see Section l2~2|) . and the turnover at g = 24.5 mag 
is due to the (much) shorter exposure time com- 
pared to the CFHTLS Deep. 

We conclude from analysis of the CFHT data 
that all the characteristics discussed in this section 
(namely a position that allows observations from 
different sites, a large number of suitable tip-tilt 
stars for laser AO, a generally low and fairly uni- 
form Galactic extinction, coupled with the small 
number of very bright stars) make this particular 
field an excellent choice for future deep extragalac- 
tic AO observations. 

5. Real-world benefits of undertaking ob- 
servations in AO-optimized fields 

How do the properties of the fields identified in 
the previous sections compare with those of exist- 
ing extragalactic deep fields? Figure O presents a 
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Fig. 7. — Upper panel: CFHT MegaCam g' and z' band mosaic of the 1° x 1° field centered on RA = 7 h 24 m , 
DEC = — 1°27'15". Note that there are only a few bright stars in this field, with the brightest one (in the 
middle of the field) at V=9.9 mag. Lower panel: Galaxy number per 0.5 mag and per 1° x 1° as a function 
of g' magnitude. Solid line represents the expected values based on CFHTLS Deep data set and triangles 
correspond to the recovered galaxy number counts. 
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Fig. 8. — Top left: Stellar surface density as a function of Galactic latitude for two sets of deep fields: 55 
existing deep fields plus 1 square degree region within the CFHTLS W2 field (blue filled circles) and 67 low 
latitude pointings we have explored in detail (red filled circles, detailed properties listed in Table [TJ. Top 
right: Average E(B — V) as a function of Galactic latitude for the same two sets of fields. Bottom right: 
Standard deviation of E(B — V) for the data in the other panels (this panel quantifies the patchiness of the 
extinction). Our chosen field (AODF) and five other fields from TablcQ]with extreme properties are labelled. 
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comparison of the average stellar density, average 
extinction, and standard deviation of the extinc- 
tion coefficient for the set of 67 55' x 55' fields 
from Table [T] to those in 55 existing deep fieldfl 
Quantities are shown as functions of Galactic lat- 
itude, with our candidate fields shown in red, and 
existing deep fields in blue. Our preferred field is 
labelled 'F' in each panel. The figure illustrates 
that our proposed fields typically have over ten 
times the stellar density of existing deep fields cou- 
pled with extinction values and extinction varia- 
tions across the field at the high end of those in 
existing deep fieldfd. How does this translate into 
practical performance benefits for undertaking AO 
observations? Are the existing fields already good 
enough? In order to investigate these questions, 
we will define a fairly generic figure of merit for 
AO observations in Section 15.11 , and compare the 
distribution of this figure of merit in our proposed 
fields to the corresponding distribution in a typical 
existing deep field (Section I5.2[) . 

For the sake of concreteness, much of the fol- 
lowing analysis will be undertaken in the context 
of the predicted performance of the soon-to-be- 
commissioned Gemini M ulti-Conjugate Adaptiv e 



Optics System (GeMS, lEllerbroek et al.l |2003) 



GeMS is intended to feed NIR instruments with 
a high Strehl ratio beam at relatively short wave- 
lengths (Strehl ratios up to ~ 50% in J— band), 
and in particular to fe ed the FLAMINGOS- 2 
NIR MOS spectrograph dEikenberrv et al.l l2004h 



and t he GSAOI imaging camera ([McGregor et al 



2004) . The benefits of undertaking observations in 



the AODF are fairly obvious for science programs 
which use imaging cameras or resolved integral 
field spectrographs, but the AODF will also be of 
considerable interest for programs of NIR MOS 
spectroscopy. Some spatially resolved kinematical 
and chemical composition information can be re- 
covered with narrow slits if these span individual 
objects that are not obliterated by seeing. Fur- 



4 From the list compiled by J. Brinchmann, see http://www. 
strw. leidenuniv. nl/~jarle/Surveys/DeepFields/ index . 



Ihtml 

5 Note that in Figure[8]we have also included a region within 
the Canada-France Hawaii Survey Legacy Survey Wide-2 
Field (CFHTLS W2). This is not technically a 'deep field', 
but CFHTLS W2 is worthy of inclusion because it is the 
most commonly-observed extragalactic field that lies in the 
general vicinity of our preferred field (for detailed descrip- 
tions and comparisons with this field, see Section l5,2j l, 



thermore, AO-assisted MOS spectrographs will be 
able to effectively use narrow slits, which minimize 
background contamination. 

We have chosen to focus our analysis on GeMS 
because it is the most advanced AO system likely 
to be available on 8m-class telescopes for the 
foreseeable future, and because Multi-Conjugate 
Adaptive Optics (MCAO) is virtually certain to 
be an important operational mode for future 30m- 
class Extremely Large Telescopes. We will not de- 
scribe the fu ndamentals of MCAO here, and refer 
the reader to iRigaut et al. ( 2000l ) for an explana- 
tion of the general principles. For our purposes it 
suffices to note that MCAO's purpose is to deliver 
high image quality over a wider area than conven- 
tional adaptive optics systems, and it does so by 
using a constellation of laser guide stars beacons 
and several natural guide stars to determine the 
shape of the wavefront, which is then corrected by 
multiple deformable mirrors. (The natural guide 
stars are essential, because they are needed in or- 
der to establish tip-tilt corrections). In the case 
of the GeMS MCAO system, five artificial guide 
star beacons and three natural guide stars are 
used. Another important point that should be 
emphasized here is that the precision of the tip- 
tilt corrections depends on distance from the guide 
stars, so the geometric arrangement of the natu- 
ral guide stars plays an important role in estab - 
lishing AO performance (jFlicker fc Rigautll200ll ). 
We assume nominal performance of GeMS, and 
emphasize that both the areal coverage and PSF 
stability expected from MCAO are substantially 
larger than in the case with conventional (single- 
laser, single wavefront sensor) AO. 

5.1. Figure of merit 

Any number of figures of merit can be devised 
for inter-comparing the performance of various AO 
systems, but in this paper we propose to use a 
figure of merit that captures the basic idea that 
real-world AO performance generally depends not 
only on image quality, but also on the variation 
of that image quality over the field of view. We 
therefore adopt the following figure of merit: 



F 



1 



T 0.25 



x (1 - 



(2) 



where (S) is the average Strehl ratio achieved in 
the field of view, and as is the RMS variation in 
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Strehl over the field of view. In the present paper 
our purpose is to understand the impact of tip- 
tilt stars, so we will be calculating Strehl ratios 
using simulation software that computes the dis- 
tortion in the wavefront due to anisokinetism and 
assumes perfect correction for other aberrations in 
the wavefront. Our procedure for doing this will 
be described in the next section. 

The distribution of F across the sky character- 
izes the performance of an AO system. The spe- 
cific values of the exponents in our definition of 
F are chosen to weight the peak Strehl at the ex- 
pense of some uniformity in the value of the Strehl 
ratio over the field of view. However, uniformity 
in Strehl is not completely de-emphasized, and a 
guide star configuration resulting in a generally 
high but strongly variable Strehl across the field 
of view will have a significantly lower value of F. 
None of main conclusions of this paper are strongly 
dependent on the specific values of the exponents 
used in Equation [2] as shown in Appendix [B] 

5.2. MCAO observations in existing and 
proposed fields 

To investigate limitations in MCAO perfor- 
mance in various fields as a function of natural 
reference star magnitudes and configurations, we 
used the Gemini MCAO simulator (F. Rigaut, pri- 
vate communication) to compute the distribution 
of the figure of merit F for ~ 1000 uniformly dis- 
tributed pointings within several fields. As noted 
earlier, we wish to study the errors introduced 
into the corrected wavefront by a paucity of tip- 
tilt stars, so the Strehl ratio used in our analy- 
sis isolates the RMS contribution to the distorted 
wavefront introduced by anisokinetism. In other 
words, the simulation assumes that all other con- 
tributions to the wavefront degradation are negli- 
gible, so if the tip-tilt correction were perfect, the 
Strehl ratio would be unity. In practise of course 
errors other than anisokinetism will contribute to 
the wavcfroni@. However, the point is that our 
analysis allows us to study the best performance 
possible from the AO system, limited only by the 
number of natural tip-tilt reference stars. 



6 The interested reader is directed to Table 2 of this web 
page for a census of other contributions to the MCAO 
wavefront: 

http: / /www. gemini . edu/sciops/ instrument s/mcao? 



Each simulated pointing was 80" x 80" in size 
(appropriate to Gemini's GeMS). We analyzed 
performance in the AODF, and for comparison 
with performance in un-optimized deep fields, we 
also examined ~ 1000 uniformly distributed posi- 
tions in the COMBO 17 Field 2 deep field and in 
the 1 square degree region of the CFHTLS wide 
fie ld (W2). COMB O 17 Field 2 (labeled as A 901 
in IWolf et all 120031 ) was chosen because it is the 
existing deep field with the highest star counts. 
We chose to include (a part of) the CFHTLS W2 
field in this comparison because (a) it is located 
fairly close to the AODF, and (b) its star counts 
surface density is high. However, large portions 
of this 7° x 7° field are 'contaminated' with very 
bright starfl In order to avoid those regions we 
performed our analysis in the 1° x 1° subfield 
within CFHTLS W2 centered on RA = 8 h 42 m , 
DEC = — 1°15', which has few very bright stars 
but an abundance of stars in the range suitable for 
tip-tilt correction. 

To determine the optimal set of guide stars 
in each pointing, we looked for the best set of 
three stars with magnitudes in the range 13 < 
R < 16.5 mag with distances between 40" and 
60" from the field center (these distances are set 
by the patrol fields of pick-off mirrors in GeMS). 
As noted by iFlicker fc Rigautl (|200lh . the ideal 
geometry for these stars is an equilateral trian- 
gle, so we searched for three stars lying in the 
40" — 60" annulus whose interior angles were 
within 60° ± 20° from each other. If we failed 
to find suitable stars we then relaxed our criterion 
that the guide stars approximate equilateral trian- 
gles and simply searched for three stars defining a 
triangle with any set of side lengths. Where multi- 
ple triangles existed we retained the one that gave 
the best value of F. When three stars in any ge- 
ometry could not be found we calculated MCAO 
performance with available stars, either singly or 
in pairs. In cases where no stars in the suitable 
magnitude range were found in the vicinity of our 
pointing, the value of F was set to zero. 

The results from our investigation are shown 
in Figure [9j which compares the distribution and 
cumulative distribution functions of F for the 
AODF, shown in blue, with the corresponding dis- 



7 See http : //legacy . astro . utoronto . ca/Fields/images/ 
w2 . html I 
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Fig. 9. — Cumulative distribution (left) and distribution (right) of our Strehl-related figure of merit for 
55' x 55' AODF (blue), COMBO 17 Field 2 (red), and the CFHTLS W2 subfield (green). Data points are 
binned according to their figure of merit value in 3 unit wide bins. The left-side vertical axis gives the 
percentage of field area with the figure of merit value equal or higher than the corresponding bin's lower 
limit. In the right-side panel ordinate represents the number of pointings contained in each bin. Both panels 
correspond to a guide star magnitude limit of 16.5 mag in i?-band (foreseen bright time limit for the Gemini 
GeMS MCAO system). See text for details. 



tributions for the COMBO 17 Field 2 and a 1 
square degree region of the CFHTLS W2 field, 
shown in red and green, respectively. The two 
panels of Figure [9] correspond to a limiting natu- 
ral guide star magnitude of 16.5 mag in R-band 
(corresponding to the anticipated limit for GeMS 
bright time observing). Although the limiting 
guide-star magnitude will vary with lunation, the 
only difference in the case where the faintest natu- 
ral guide stars is R — 18.5 mag (the ideal dark sky 
performance of GeMS) will be a higher number of 
pointings with low F. 

The enormous benefits that emerge from ob- 
serving in an AO-optimized field are obvious from 
even a cursory inspection of Figure El but it is 
worth discussing the figure in some detail. In 
analyzing this figure, we noted that even in our 
adaptive-optics optimized deep field (ID 8328, see 
Table [I]), for bright time observing only about 
~ 1% of the pointings in the field fall within tri- 
angles of guide stars that are approximately equi- 
lateral (as defined above). However, if we allow 
any configuration of three suitably bright guide 
stars, sky coverage in the AODF is 55%. The 



^ F < 20 range in this case covers the MCAO 
performance when less than three natural guide 
stars are available. If we assume dark time ob- 
serving, the sky coverage is ~ 92% but the high- 
est figure of merit still corresponds to the con- 
figurations of brighter (R ^ 16.5) stars. In con- 
trast to this, none of ~ 1000 uniformly distributed 
pointings in the COMBO 17 Field 2 and only one 
pointing in the CFHTLS W2 subfield fell within 
equilateral triangles of natural guide for our sim- 
ulated bright time conditions, and this fraction 
rises to only 0.2% (COMBO 17 Field 2) and 1.2% 
(CFHTLS W2) at dark time. If the MCAO re- 
quirement is relaxed to allow any configuration of 
three NGS, COMBO 17 Field 2 coverage is 0.4% 
and 3.8% for bright and dark time observing con- 
ditions, respectively. Three AO-suitable stars with 
R ^ 16.5 forming a triangle were found around 
1.8% of the CFHTLS W2 pointings. If the mag- 
nitude limit is lowered to R < 18.5, CFHTLS W2 
coverage for sets of three stars arranged in any 
type of triangle reaches 21%. For both COMBO 17 
and CFHTLS W2, the most common figure of 
merit values that emerge (0 < F < 20) correspond 



17 



1 

10" 2 
lO" 4 



npj iiiiii^ ■ ill iii^ iiiiii^ i iiing^^mi 1 1 mm J p 



10 



i| iiiiiiJ iiiiiii| iiiiiiJ iiiiiiiI ^ i ^ mi 



10" 



io -2 L. 



IO" 4 ;. 




O s 



AODF 



10" 2 ^___^_ # __ j 



io -4 - 




i| iiiiiij iiiiiii| iiiiiiJ iiiiiii| ^ i mj 1 1 il l ! £ 



$CFHTLS W2l 

o 




IO" 1 _ 



10" 



10"' _ 



10 "iJ j J ^ I A 1 

10" 6 10" 4 10" 2 



AODF 



I I III LB 



° COMBO 17 
i 1 1 m ill m i m ill i 1 1 ii ui 

° a >*b(g > 1 



o CFHTLS W2 



10" l 
1 IO" 3 



Q3 
M 

Qi 



:82 

1425 

375 
325 
275 
225 
175 

125 

78.7 

44.6 
16.7 



10" 



10" 1 



1 



FoM 



Fig. 10. — Strehl ratio as a function of its RMS variance for pointings accessible to AO in AODF, COMBO 17 
Field 2, and CFHTLS W2 1 square degree subficld. Left-side panels display full range of values for Strehl 
ratio and its variance. The dashed line corresponds to MCAO 'failure' for pointings where only one or two 
guiding stars are available. The right-side panels show the range of Strehl ratios and related variances for 
pointings where MCAO is efficient (i.e., where at least one set of three AO-friendly stars forming a triangle 
is available). In each panel points are color-coded according to figure of merit value (equation 2). Color 
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distribution across the AODF. All panels correspond to a guide star magnitude limit of 16.5 mag in i?-band 
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to MCAO performance with only one or two stars 
natural tip-tilt stars. Thus the low sky coverage 
in the COMBO 17 Field 2 and CFHTLS W2 sub- 
field results in the large spike in the first bin in 
Figure |H1 with only a few (out of 1000) positions 
with high F producing a steep decline at F f=s 20. 
In contrast to this, the median value of F in the 
AO-optimized field is a factor of three higher than 
in COMBO 17 Field 2 and CFHTLS W2, and a 
significant tail of F extends out to F ~ 500 (a 
factor of three better than the best F obtained in 
COMBO 17 Field 2 and CFHTLS W2 subfield). 

Figure [9] shows the superiority of the AODF 
over both COMBO 17 Field 2 and CFHTLS W2 in 
terms of the figure of merit. The relation between 
the figure of merit values and the quantities that 
enter Equation [2] (the average Strehl ratio and its 
RMS variation over the field of view) is explored 
in Figure [T0l where points are color-coded accord- 
ing to their figure of merit values. (The reader 
is once again reminded that the Strehl ratios in 
our calculations take into account only the tip- 
tilt correction). Dashed line in the left-hand pan- 
els denotes threshold value of the average Strehl 
ratio at which MCAO fails, i.e., when there are 
only one or two NGSs available. While ~ 70% 
of the AODF pointings with available NGS have 
Strehl ratios above the MCAO threshold, for the 
COMBO 17 Field 2 and the CFHTLS W2 that 
fraction is 3.8% and 8%, respectively. The right- 
hand panels show the distribution of Strehl ra- 
tios and corresponding variances for MCAO per- 
formance when the required geometry of guiding 
stars is available. Figure ITUl confirms that the fig- 
ure of merit we have defined does rather nicely 
map onto fields with the combination of high aver- 
age Strehl ratios and low RMS variation in Strehl 
over the 80" x 80" field of view. The third and 
fourth quartile of the figure of merit distribution 
for the AODF (orange and red points in Fig. fTO)) , 
that enclose Strehl ratios of 90% < (S) < 97% 
and RMS variations of 1% < as < 6%, contains 
only 3% of the CFHTLS W2 pointings available 
for AO and 1.5% of the corresponding pointings 
within the COMBO 17 Field 2. This hugely bet- 
ter AO performance clearly illustrates the bene- 
fits of undertaking extragalactic AO observations 
in fields optimized for adaptive optics. 



5.3. The AODF in the era of upcoming 
methods and instruments 

Although the MCAO has been emphasized in 
the present paper, other AO methods are in de- 
velopment especially for use with the 30m-class 
telescopes. These methods include multi-object 
adaptive optics (MOAO) and systems with faint 
infrared (IR) tip-tilt sensors whose images are 
"sharpened" by the AO system. Future AO sys- 
tems are expected to be less sensitive to the 
(bright) NGS surface density. 

MOAO (jAssemat et all 120071 ) is a technique 
that allows simultaneous AO corrections for sev- 
eral small IFU target fields (typically 2 — 4" in 
diameter, sufficient to map velocity fields of large 
spiral galaxies at z > 1) within a wider field of 
view (FOV^ 5 — 10', driven by the surface den- 
sity of line emitting galaxies at z > 1). Each 
IFU target field is corrected by a separate de- 
formable mirror (DM) that provides AO correc- 
tion along a given line-of-sight only (in contract 
with MCAO systems that provide this correction 
across the whole FOV) . Multiple NGS are used for 
tomo graphic wavefront sensing (jRagazzoni et al^ 



2000), i.e., to probe the 3D phase perturbations in 
the atmosphere above the telescope primary aper- 
ture. A real-time control system then slices multi- 
ple columns through the mapped turbulent volume 
in the directions of all targets and applies a cor- 
recting signal to multiple independent DMs. The 
critical difference between MOAO and MCAO is 
that the former is an open loop system where the 
wavefront sensors do not get any feedback from 
DMs. In other words, the wavefront error is mea- 
sured and corrected only once, and the accurate 
system calibration is essential. 

While an advantage of open-loop MOAO is 
that it does not need a guide star for each tar- 
get, the limiting magnitude of the NGS is set 
by the requirement for low system error level. 
Thus even MOAO designed for the 30m-class tele- 
scopes will still need guide stars brighter than 
R = 17 mag. For example, an appropriate con- 
figuration for mapping turbulent atmospheric lay- 
ers using the EAGLE MOAO system that is be- 
ing developed for the E-ELT involv es 5 NGS with 



R < 17 mag in the 7f 3 patrol FOV (jRousset et al 
2010). Based on the average star counts surface 
density, only one out of 55 existing deep fields 
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mentioned earlier (COMBO 17) lies above (and 
very close to) this threshold. On the other hand, 
the number density of available NGS in the AODF 
is six times higher than EAGLE MOAO require- 
ments. 

As was noted earlier, as the primary mirror of 
a telescope becomes larger, its sensitivity allows 
fainter stars to be used for wavefront sensing. AO 
systems on 30m-class telescope will give satisfy- 
ing performance with a NGS magnitude limit of 
R f« 21 mag (though the limiting magnitude will 
be lower in MOAO mode, as noted above). We 
performed the same analysis as in previous sec- 
tion on the COMBO-17 field using GeMS MCAO 
simulator and looked for stars arranged into trian- 
gles within magnitude range 13 < R < 21, based 
on the US NO-B catalog (with limiting magnitude 
of R = 22, iMonet et al1l2003l ). Although the sky 
coverage in this case is ~ 47%, only ~ 15% of 
these sets of three stars will provide Strehl ratio 
greater than 0.5. On the other hand, the sky cov- 
erage for high Strehl ratio values (i.e., guide stars 
arranged in equilateral or isosceles triangles) using 
MCAO system in the AODF is 100%. 

Another advanced type of AO planned for 30m- 
class telescopes will utilize IR tip-tilt wavefront 
sensing. A major advantage of this approach is 
increased guide star density, since faint IR NGS 
images are sharpened by the AO system. For 
example, for the TMT NFIRAOS the probabil- 
ity of finding at least 1 tip-tilt star brighter than 
J = 21 mag is 95% at high galactic latitude. How- 
ever, at least three NGS are still required to de- 
tect the effects of tilt anisoplanatism; the use of 
only one off-axis tip-tilt star would give blurred 
time-averaged images of the science objects. This 
condition lowers the sky coverage tha t NFIRAOS 



will a chieve to 50% at high latitudes (jWang et al 



2008). Even with three guide stars one expects to 
get a continuum of performance. Fainter stars will 
force the system to run slower which in turn leaves 
larger tip, tilt and focus errors. Thus, although 
MCAO system on a 30m-class telescope will oper- 
ate over much of the sky, an insufficient number 
of bright guide stars will impair imaging perfor- 
mance, as diffraction-limited cores will be blurred 
out by these tip/tilt /focus errors and the variation 
across the FOV will be increased. (IFU work will 
be less affected, because the ensquared energy loss 
in a spaxel a few times larger than the diffraction 



limit will be lower) . This type of systems will give 
the best results when used on a field densely popu- 
lated with bright NGS, where the AO is not push- 
ing the boundaries of the control system. Finally, 
we note that AO using natural guiding stars only 
(of magnitude ~ 12 and within ~ 15"of the science 
object ) is capable of achieving higher Strehls than 
MCAO, is easier to do, and removes all the compli- 
cations of changing plate scales, cone effect, laser 
elongation, etc. A major benefit of doing AO in 
the proposed field is that it will allow much of the 
of the 30m-class telescope science to be done in 
NGS mode. For example, if an AO system, simi- 
lar to EAGLE MOAO, that uses only bright NGS 
were employed in the AODF, sky coverage would 
be > 75%. 

We conclude this section by noting that fields 
optimized for ground-based adaptive optics with 
30m-class telescopes will not be made obsolete 
by upcoming space missions, such as JWST and 
Euclid/WFIRST (both of which are 5-10 years 
away in any case). Euclid/WFIRST will likely 
only operate in slitless spectroscopy mode, while 
JWST will be equipped with a micro-shutter ar- 
ray for simultaneous spectroscopy of ~ 100 sources 
and with an IFU for 3D spectros copy, spanning 
the w avelength range of 1 — 5 fj,m (jGardner et al 



2006). However, the main point is that the future 
ground-based observations will still be undertaken 
at spatial resolutions a factor of 5 and 15 times 
higher than the angular resolution of JWST and 
Euclid/WFIRST, respectively. 

6. Summary and Conclusions 

We have combed through stellar density and 
extinction maps to identify 67 low Galactic lati- 
tude fields with high star density, remarkably low 
extinction, relatively large area (1 square degree, 
to mitigate the effects of cosmic variance) , and an 
AO-friendly stellar mix of many R = 13 — 16.5 mag 
but few < 8 mag stars. These fields allow highly 
efficient adaptive optics to be undertaken in low 
extinction extragalactic fields with minimal sat- 
uration and scattering. A comparison of these 
fields with existing deep fields reveals that while 
the number of guiding stars per square arcminute 
is on average 15 times higher in the AO-friendly 
fields, the mean level of extinction is comparable 
to the more extincted existing deep fields. 
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By augmenting our analysis with some practical 
considerations (such as the desire for an equatorial 
field accessible from both hemispheres), we iden- 
tify a single one square degree field (which we des- 
ignate the Adaptive Optics Deep Field, or AODF) 
as being particularly promising for extragalactic 
AO work. This field is centered at RA: 7 h 24 m 3 s , 
Dec: — 1°27'15". Analysis of galaxy counts in this 
field based on short observations of this AODF in 
g' and z' bands (using MegaCam on CFHT) con- 
firm both the absence of extinction and the abun- 
dance of suitable tip-tilt stars the AODF. In fact, 
galaxy counts in this field closely follow the counts 
found in the CFHTLS Deep data set. 

Simulations were undertaken to estimate the 
practical performance benefits of undertaking AO 
observations in the AODF. Our analysis shows 
enormous advantages emerge from undertaking 
AO observations in optimized fields such as the 
one described here. For example, for geometries 
of natural guide stars which produce spatially sta- 
ble high Strehl ratio PSFs, dark time sky coverage 
in the AODF is essentially 100% using the Gemini 
MCAO system, which is a factor of over 50 times 
higher than for most existing deep fields. 
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Appendices 



A. Properties of selected fields 

E(B— V) Star counts surface density 





0.50 . 



rZ_ i.U jj/muT 
(275.9, -41.1) Equatorial (275.9, -41.1) Equatorial 

Fig. 11. — Extinction and star surface density maps of the 6'871 x 6'871 cells within the 55' x 55' held with 
the highest mean star counts surface density of (E sc ) = 1.3 arcmin -2 (i.e., with 11,000 star-forming galaxies 
potentially observable with AO if the average number surface density of these objects is 3 arcmin -2 ). The 
held is centered at a = 18 h 23 m 39?53, 5 = — 41°6'44'.'3. Solid line grid corresponds to the celestial coordinate 
system with RA[°] = 360° - a[°]. 

Figure @] demonstrates how the set of 55' x 55' fields that are most suitable for AO observations were 
selected. Amongst this set we chose our preferred field (the 'AODF') partly on the basis of accessibility 
to both hemispheres. Other fields may be equally (or even more) suitable if this criterion is relaxed. A 
field of obvious interest is the one with the highest star count surface density. Table [1] shows that this field 
is located at a = 18 h 21 m 57?52, S = -43°14'2'.'775. It is flagged with an open circle in Figures H and M 
and presented in more detail in Figure 1 1 1 1 where star count surface density and extinction maps with finer 
sampling (6'871 x 6f 871) are used to bring out the features of individual cells within the field. The number 
of stars in this field is ~ 5000, corresponding to more than one star per arcminute squared. 

Another property worth optimizing for is dust extinction, so it is interesting to look for fields with 
extraordinarily low extinction in Table 1. The field with the lowest reddening coefficient E(B — V) is 
centered at a = 18 h 30 m 14!96, S = -47°35'20'.'76 and also flagged in Figures H and [8] A close-up view of 
its stellar surface density and reddening coefficient distribution is given in Figure [T^l Although this field 
features lower extinction then some of already explored deep fields, the number of its potential AO guiding 
stars is almost an order of magnitude higher then in the existing deep fields labeled in Figure [8] 

Finally, it is interesting to consider fields with highly homogenous extinction. The third flagged field 
in Figures 0] and [8] is the one with the lowest value of standard deviation for E(B — V) from Table [1] at 
a = 17 h 40 m 30!18, 5 = -61°25'54'.'54. Tthe properties of this field's 6.'871 x 6f871 cells are given in Figure[Lj 
Despite not having the highest number of AO suitable stars (> 2350 vs. > 4600) or the lowest dust extinction 
(its (E(B — V)) is ~ 35% higher then in the field with minimum extinction), this field may be interesting 
for certain studies in which reddening homogeneity is more important than other factors. 

The higher order statistical moments of star counts within sub-fields, along with proximity to other 
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Fig. 12. — Extinction and star surface density maps of the 6f871 x 6.'871 cells within the 55' x 55' field 
with the lowest reddening coefficient ((E(B — V)} — 0.053 mag) but still containing more then 2500 possible 
guiding stars (~ 5900 observable star- forming galaxies with number surface density of 3 arcmin -2 ). The field 
is centered at a = 18 h 30 m 14 i ;96, 5 = — 47°35'20'.'76. Solid line grid corresponds to the celestial coordinate 
system with RA[°] = 360° - a[°}. 
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Fig. 13. — Extinction and star surface density maps of the 6f871 x 6'871 cells within the 55' x 55' field with 
the lowest standard deviation of E(B — V) that translates into the smallest variation in dust content across 
the field. There are > 2350 13 — 16 mag stars suitable for AO guiding in the field corresponding to ~ 5600 
observable star- forming galaxies with the number surface density of 3 arcmin -2 . The field is centered at 
a = 17 h 40 m 30?18, S = — 61°25'54'.'54. Solid line grid corresponds to the celestial coordinate system with 
RA[°] = 360° -a[ }. 
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suitable areas, may be other factors worth considering when choosing fields. Whether the the higher order 
moments really matter depends on the specific science objectives of the observations. The skewness, 73, in 
the star count surface density distribution might be worthwhile to consider in cases where one wishes to 
optimize for having a smaller number of fields with many NGS. For example, 73 < corresponds to the mass 
of distribution shifted towards higher values. For such fields there are many patches of very high star density. 
On the other hand, kurtosis in the E(B — V) distribtuion, 74, could be important if one wishes to optimize 
a field for photo-z consistency. Fields with 74 < have a less peaked distribution of E(B — V), i.e., more 
uniform extinction. The area within dashed lines in Figure |4] contains the fields from three regions in the 
sky. When identifying the best AO field in each region (coloured arrows in Figures |4] and [5]), we have taken 
into account all four moments of both star counts surface density and extinction distributions. In Figure [5] 
we present three moments - standard deviation, skewness, and kurtosis - for all 67 fields from Table Q] as 
functions of the mean star count surface density (the most important factor for identifying AO-friendly 
fields). 

B. The form of the figure of merit 




30 40 100 300 500 1000 



i/« X (i-<s>r) 

Fig. 14. — Cumulative distribution of our Strehl- related figure of merit for 55' x 55' AODF (black), 
COMBO 17 (red) field, and the 1° x 1° subfield of CFHT W2 (green). Data points are binned accord- 
ing to their figure of merit value in 3 (units) wide bins. Vertical axis gives the percentage of field area 
with the figure of merit value equal or higher than the corresponding bin's lower limit. In each panel the 
figure of merit is defined by different combination of the two exponents (Eq. IB1[) . All panels correspond to a 
guide star magnitude limit of 16.5 mag in i?-band (foreseen bright time limit for the Gemini GeMS MCAO 
system). 

Our proposed figure of merit, given in Equation[2l is useful for characterizing AO system performance, but 
there is considerable flexibility in choosing the values of the exponents in this equation. We chose exponents 
which strike a balance between emphasizing the importance of peak Strehl ratio in a field and emphasizing 
the uniformity of the Strehl ratio throughout the field. Other, equally valid, choices of these exponents could 
be made that strike a different balance. In order to investigate how different combinations of exponents in 
Equation IB 11 might influence our general conclusions, we have defined a more generic form of the figure of 
merit: 
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Fig. 15.— Distribution of our Strehl-related figure of merit for 55' x 55' AODF (black), COMBO 17 (red) 
field, and the 1° x 1° subfield of CFHT W2 (green). Data points are binned according to their figure of merit 
value in 3 (units) wide bins. Ordinate represents the number of pointings contained in individual bin. In 
each panel the figure of merit is defined by different combination of the two exponents (Eq. IBlj) . All panels 
correspond to a guide star magnitude limit of 16.5 mag in i?-band (foreseen bright time limit for the Gemini 
GeMS MCAO system). 



F a%x(l-(S))f>> (B1) 

where exponents a and /3 can take several values: a £ {0.25,0.125}, (3 £ {1, 1.5,2}. The maximum image 
quality (i.e., average Strehl ratio) is weighted by the (3 exponent. Constraints on the range of values for a 
are set so that the high values of the figure of merit cannot coincide with extremely low values for both Strehl 
ratio and its variance. 

The resulting distributions for all six possible combinations are given in Figures [U and 1151 In this figure 
we have examined and compared three fields: AODF, COMBO 17, and a 1° x 1° subfield within CFHT W2. 
The corresponding histograms are presented in blue, red, and green, respectively. Of course the range of 
values that the figure of merit can take is seen to depend rather strongly on the choice of exponents, but the 
important thing to note is that none of the distributions show major changes in shape or relative position 
for different combinations of a and fi. Furthermore, the highest figures of merit in all six cases (for all three 
fields) correspond to the pointings with the highest average Strehl ratio ((S) > 90%) and relatively low 
variation of Strehl across the 80" x 80" field of view (erg < 2%). Since our analysis relies only on relative 
comparisons between different fields, our overall conclusions seem quite robust to the specific choices of the 
exponents. 
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Table 1 

Properties of the 55' x 55' fields identified in Figures [3] & [4] 
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a pixel number based on iV gidc = 64 HEALPix partition of the sphere 

the third moment of distribution (skewness) for star counts surface density 
the fourth moment of distribution (kurtosis) for star counts surface density 
the third moment of distribution (skewness) for extinction 

e the fourth moment of distribution (kurtosis) for extinction 
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